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We determine masses and decay constants of heavy-heavy and heavy-charm pseudoscalar mesons 
as a function of heavy quark mass using a fully relativistic formalism known as Highly Improved 
Staggered Quarks for the heavy quark. We are able to cover the region from the charm quark 
mass to the bottom quark mass using MILC ensembles with lattice spacing values from 0.15 fm 
down to 0.044 fm. We obtain /b, = 0.427(6) GeV; ms, = 6.285(10) GeV and = 0.667(6) 
GeV. Our value for f^^ is within a few percent of fr confirming that spin effects are surprisingly 
small for heavyonium decay constants. Our value for /s^ is significantly lower than potential model 
values being used to estimate production rates at the LHC. We discuss the changing physical heavy- 
quark mass dependence of decay constants from heavy-heavy through heavy-charm to heavy-strange 
mesons. A comparison between the three different systems confirms that the Be system behaves in 
some ways more like a heavy-light system than a heavy-heavy one. Finally we summarise current 
results on decay constants of gold-plated mesons. 



I. INTRODUCTION 

Lattice QCD calculations offer particular promise for 
B meson physics where a number of relatively simple 
weak decay processes give access to elements of the CKM 
matrix that are important for constraining the unitarity 
triangle of the Standard Model pj . The theoretical cal- 
culation of the appropriate weak matrix elements must 
be done with percent accuracy for stringent constraints, 
making optimal use of the experimental results. This has 
not yet been achieved, despite the enormous success of 
lattice QCD over the last five years and its acceptance 
as a precision tool for QCD physics [2]. Work is ongoing 
on several different approaches. Here we continue discus- 
sion of an alternative method for B meson physics that 
may offer a faster route to high accuracy for some quanti- 
ties than other methods currently in use. Following work 
on accurate b and c quark masses |3] and heavy-strange 
decay constants [4], we show results for masses and de- 
cay constants of Be and rji, mesons and map out their 
heavy-quark mass dependence. As well as showing that 
high accuracy can be achieved, these results provide an 
interesting comparison of how heavy-charm mesons sit 
between heavy-heavy and heavy-light. 

The calculations use a discretisation of the quark La- 
grangian onto the lattice known as the Highly Improved 
Staggered Quark (HISQ) action This has the advan- 
tages of being numerically very fast along with having 
small discretisation errors and enough chiral symmetry 
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(a PCAC relation) that the weak current that causes 
charged pseudoscalar mesons to decay leptonically is ab- 
solutely normalised. This action readily gives tt and K 
meson decay constants with errors below 1% on giuon 
field configurations that include the full effect of u, d and 
s quarks in the sea |BJ |7] . Results from multiple values 
of the lattice spacing and multiple sea u/d quark masses 
allow extrapolation to the real world with physical u/d 
quark masses at zero lattice spacing. 

The HISQ action gives similarly accurate results for 
mesons containing c quarks [6] , significantly improving on 
previous methods that use a nonrelativistic effective the- 
ory such as the Fermilab action or NRQCD (Hj. The 
key advantages are clear: the HISQ action has no errors 
from missing higher order terms in the effective theory 
or from the renormalisation of the decay constant [T]. 
The price to be paid is that of the discretisation errors. 
These errors are much larger for c quarks than for u/d 
and s, since their size is now set by mcO, rather than 
Aqcdo-- They can be well controlled, however, using the 
HISQ action on gluon configurations with a wide range 
of lattice spacing values down to 0.045 fm where TOcO = 
0.2 [TO]. Discretisation errors are in fact the only issue 
for the Dg meson, for which particularly accurate results 
can be obtained. This meson has no valence light quarks 
and the dependence of both its mass and decay constant 
on the u/d quark masses is seen to be very small |10) . 
meaning that uncertainties from the chiral extrapolation 
are not significant. 

It is less clear what to do for b quarks because they are 
so much heavier. To achieve m^a < 1 we need a lattice 
spacing, a < 0.04fm. Using NRQCD or the Fermilab for- 
malism we can readily handle b quarks on much coarser 
lattices, with a ss 0.1 fm, but must then take a substan- 
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tial error (currently 4% for NRQCD [TP) from matching 
the weak annihilation current to full QCD perturbatively. 
Work is underway to reduce this error [H] . It should also 
be emphasised that this matching error is not present in 
ratios of decay constants, for example fss/ fs which is 
known to 2% from NRQCD [11] ■ 

Here we show what accuracy is possible using the HISQ 
action for b quarks. We use quark masses heavier than 
that of the c quark and map out the heavy quark mass 
dependence of both masses and decay constants for a va- 
riety of different pseudoscalar mesons. By using experi- 
ence from the Dg [10] and concentrating on mesons that 
do not contain valence light quarks we do not have to 
worry significantly about the extrapolation to the phys- 
ical u/d quark mass limit. The key issue is that of dis- 
cretisation errors, as for f^^ , and we therefore work with 
the same large range of lattice spacing values from 0.1 5fm 
to 0.044fm, so that we can account fully for the a depen- 
dence. It is important to separate discretisation effects 
from physical dependence on the heavy quark mass since 
we do also have to extrapolate to the physical b mass from 
the quark masses that we are able to reach on these lat- 
tices. We are only able to obtain results directly at close 
to the physical b mass on the finest, 0.044fm lattice. 

We have already demonstrated how well this method 
works in determining the decay constant of the Bs me- 
son [4], one of the key quantities of interest for CKM 
studies. Mapping out the Bg decay constant as a func- 
tion of heavy quark mass showed that the decay constant 
peaks around the Dg and then falls slowly. We found that 
fss/ — 0.906(14), the first significant demonstration 
that this ratio is less than 1. 

Here we extend this work to map out results for the 
decay constants of the 77;, and B^. mesons, along with the 
Be meson mass. The Be meson mass is known exper- 
imentally but its leptonic decay rate has not yet been 
measured and so we provide the first prediction of that 
in full lattice QCD. The masses and decay constants also 
reveal information about the nature of these mesons that 
can provide useful input to model calculations. For ex- 
ample, does the Be meson look more like a heavy-heavy 
meson or a heavy-light meson? It is important to em- 
phasise that both the results determined at the b quark 
mass and the dependence on the heavy quark mass (and 
on any light quark masses) have physical meaning: the 
former can be tested against experiment but the latter 
can provide stringent tests of models and comparison be- 
tween lattice QCD calculations. 

The layout of the paper is as follows: secti on [H] de- 
scribes the lattice calculation and then section |III| gives 
results for heavy-heavy and heavy-charm mesons in turn. 
We compare the Be meson mass to experiment and pre- 
dict its decay constant as well as comparing the be- 
haviour of heavy-charm mesons to that of heavyonium 
and heavy-strange mesons. Section |Vj gives our conclu- 
sions, looking forward to what will be possible for b quark 
physics on even finer lattices in future. 
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Nco,if X Nt 


1 


2.152(5) 0.0097 


0.0484 


16 


48 


631 X 2 


2 


2.618(3) 0.01 


0.05 


20 


64 


595 X 2 


3 


3.699(3) 0.0062 


0.031 


28 


96 


566 X 4 


4 


5.296(7) 0.0036 


0.018 


48 


144 


201 X 2 


5 


7.115(20) 0.0028 


0.014 


64 


192 


208 X 2 



TABLE I; Ensembles (sets) of MILC configurations used for 
this analysis. The sea asqtad quark masses m""'' (l — u/d) 
and rris'"' are given in the MILC convention where uq is the 
plaquette tadpole parameter. The lattice spacing values in 
units of ri after 'smoothing' are given in the second col- 
umn [T3]. Set 1 is 'very coarse'; set 2, 'coarse'; set 3, 'fine'; set 
4 'superfine' and set 5 'ultrafine'. The size of the lattices is 
given by X T. The final column gives the number of config- 
urations used and the number of time sources for propagators 
per configuration. 



II. LATTICE CALCULATION 

We use ensembles of lattice gluon configurations at 5 
different, widely separated, values of the lattice spacing, 
provided by the MILC collaboration. The configurations 
include the effect of u, d and s quarks in the sea with the 
improved staggered (asqtad) formalism. Table |l] lists the 
parameters of the ensembles. The u and d masses are 
taken to be the same, and the ensembles have rriu/d/'^s 
approximately 0.2. As discussed in section |lj we expect 
sea quark mass effects to be small for the gold-plated 
mesons with no valence light quarks that we study here. 

The lattice spacing is determined on an ensemble-by- 
ensemble basis using a parameter ri that comes from fits 
to the static quark potential calculated on the lattice [13] . 
This parameter can be determined with very small sta- 
tistical/fitting errors. However, its physical value is not 
accessible to experiment and so must be determined us- 
ing other quantities, calculated on the lattice, that are. 
We have determined ri = 0.3133(23) fni using four differ- 
ent quantities ranging from the (2S-1S) splitting in the 
T system to the decay constant of the rjs (fixing fx and 
/tt from experiment) |14j . Using our value for ri and the 
MILC values for ri/a given in Table |l] we can determine 
a in fm on each ensemble or, equivalently, in GeV 
needed to convert lattice masses to physical units. It 
is important to note that the relative values of a (from 
ri/a) are determined more accurately than the absolute 
values of a (from ri). Our fits account for this to give 
two separate errors in our error budgets. 

Table |T| lists the number of configurations used from 
each ensemble and the number of time sources for the 
valence HISQ propagators per configuration. To increase 
statistics further we use a 'random wall' source for the 
quark propagators from a given time source. When quark 
propagators are combined this effectively increases the 
number of meson correlators sampled and reduces the 
statistical noise by a large factor for the case of pseu- 
doscalar mesons. We also take a random starting point 
for our time sources for the very coarse, coarse and fine 
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ensembles. 

We use many different masses for the HISQ valence 
quarks varying from masses close to that of the s quark 
to much heavier values for c quarks and for quarks with 
masses between c and b. On all sets the largest valence 
quark mass in lattice units that we use is m/jO = 0.85. 
These propagators are combined to make goldstone pseu- 
doscalar meson correlators at zero momentum with all 
possible combinations of valence quark masses. We sep- 
arate them into 'heavy-heavy' correlators when both 
masses are the same and are close to charm or heav- 
ier; 'heavy-charm' when one mass is close to charm and 
the other is heavier and 'heavy-strange' when one mass 
is close to strange and the other is close to charm or 
heavier. 

The meson correlation function is averaged over time 
sources on a single configuration so that any correlations 
between the time sources are removed. Autocorrelations 
between results on successive configurations in an ensem- 
ble were visible by binning only on the finest lattices. We 
therefore bin the correlators on superfine and ultrafine 
lattices by a factor of two. 

The meson correlators are fit as a function of the time 
separation between source and sink, t, to the form: 
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(1) 



for the case of equal mass quark and antiquark. i — 1 
is the ground state and larger i values denote radial or 
other excitations with the same J^^ quantum numbers. 
T is the time extent of the lattice. For the unequal mass 
case there are additional 'oscillating' terms coming from 
opposite parity states, denoted ip: 



C(i)= Va,e-^^'* + (-l)'a 



+ {t^T-t) (2) 



To fit we use a number of exponentials i, and where 
appropriate ip, in the range 2-6, loosely constraining the 
higher order exponentials by the use of Bayesian pri- 
ors [15) . As the number of exponentials increases, we 
see the value fall below 1 and the results for the fit- 
ted values and their errors for the parameters for the 
ground state i = 1 stabilise. This allows us to determine 
the ground state parameters Oi and Mi as accurately 
as possible whilst allowing the full systematic error from 
the presence of higher excitations in the correlation func- 
tion. We take the fit parameters to be the logarithm of 
the ground state masses Mi and Mi^ and the logarithms 
of the differences in mass between successive radial ex- 
citations (which are then forced to be positive). The 
Bayesian prior value for Mi is obtained from a simple 
'effective mass' in the correlator and the prior width on 
the value is taken as a factor of 1.5. The prior value for 
the mass splitting between higher excitations is taken as 
roughly 600 MeV with a width of 300 MeV. Where oscil- 
lating states appear in the fit, the prior value for Mi^ is 
taken as roughly 600 MeV above Mi with a prior width of 



300 MeV and the splitting between higher oscillating ex- 
citations is taken to be the same as for the non-oscillating 
states. The amplitudes a, and are given prior widths 
of 1.0. We apply a cut on the range of eigenvalues from 
the correlation matrix that are used in the fit of 10"^ 
or 10^'*. We also cut out very small t values from our 
fit, typically below 3 or 4, to reduce the effect of higher 
excitations. 

The amplitude, ai, from the fits in equations (Tj) 
and ([2| is directly related to the matrix element for the lo- 
cal pseudoscalar operator to create or destroy the ground- 
state pseudoscalar meson from the vacuum. Using the 
PCAC relation this can be related to the matrix element 
for the temporal axial current and thence to the decay 
constant. The PCAC relation guarantees that no renor- 
malisation of the decay constant is needed. We have: 



fp = {ma + mb). 



l2ai 
Mf' 



(3) 



for meson P. Here rua and rrih are the quark masses used 
in the lattice QCD calculation. 

fp is clearly a measure of the internal structure of a 
meson and in turn is related, for charged pseudoscalars 
such as the tt, K, D, Ds, B and Be mesons, to the ex- 
perimentally measurable leptonic branching fraction via 
a W boson: 

B{P ^ /fz 7 = ^' fpmfmp 1 ^ 

STT V nip 

(4) 

up to calculable electromagnetic corrections. Vab is the 
appropriate CKM element for quark content ab. Tp is the 
pseudoscalar meson lifetime. For neutral mesons there is 
no possibility to annihilate to a single particle via the 
temporal axial current. However, in the Standard Model 
the Bg and B are expected to annihilate to iJ-'^fJ-'' with 
a rate that is proportional to fp\V^lVtq\'^ via 4-fermion 
operators in the effective weak Hamiltonian [TS] . For the 
heavy-heavy pseudoscalar, the decay rate to two photons 
is related to its decay constant but only at leading order 
in a nonrelativistic expansion. In section [HT] we compare 
the pseudoscalar decay constant to that of its associated 
vector meson, determined directly from its decay to lep- 
tons. 

The results for masses and decay constants from fits 
in eqs. ([T]) and ([2| and using eq. ^ are in units of the 
lattice spacing, and given in this form in the tables of sec- 
tion |III| To convert to physical units, as discussed earlier, 
we determine the lattice spacing using the parameter ri. 

We then fit the results in physical units as a function 
of heavy quark mass to determine the heavy quark mass 
dependence and the physical value at the b quark mass. 
Because the bare heavy quark mass used in the lattice 
action runs with lattice spacing we need a proxy for it 
that is a physical quantity, such as a meson mass. In [4] 
we used the heavy-strange pseudoscalar mass since we 
were focussing on heavy-strange mesons. Here we choose 
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-0.1(1) 








-hl.3(7) 







form of /o 



A 



TABLE II: Estimates of shifts in MeV to be applied to the 
masses determined in lattice QCD to allow for missing elec- 
tromagnetism, c quarks in the sea and annihilation to gluons 
for the and 7)c mesons [TI^[1S]. The electromagnetic shift is 
estimated from a potential model for r]h^ rjc and Be and from 
a comparison of charged and neutral meson masses for Bs and 
Ds- The c-in-sea and gluon annihilation shifts are estimated 
from perturbation theory. The errors on the shifts are given 
in brackets. Note that for the r/s there are no shifts because 
the mass is fixed in lattice QCD |14] . 
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A{M/Mof 
A{M/Mof 

Ai av(M) \-2/9( M \b 

yl((M - M^J/Mo)' 
ayiM) -.-2/9/ • 



COOOO Cijkl 



A(M/Mo)'' 



0±10±2 1 0±4.5 

1 1 0± 2 0± 1.5 

-0.5 0±2 1 0±1.5 

1 1 0± 2 0± 1.5 

-0.5 0±2 1 0±3 

1 0± 2 0± 1.5 



TABLE III: The functional form for /o(M), the leading power 
dependence on the heavy quark mass, used in fitting the dif- 
ferent quantities described in section III using equation 
The third and fourth columns give the prior values and widths 
for the parameters A and b. In most cases b was fixed and 
then a single number is given. Likewise the sixth column 
gives the prior value and width for the djki where the sum 
was normalised so that coooo was set equal to 1. 



lattice spacing for each set of results f{M, a), we use 



the mass of the heavy- heavy pseudoscalar meson, rjfi, to 
provide the same a;-axis for all of our plots shovifing de- 
pendence on the heavy quark mass. The positions of c 
and b on these plots are then determined by the values 
of the ?7c E^nd rji, masses. 

The experimental results for the rif, and rye meson 
masses are 9.391(3) GeV and 2.981(1) GeV respec- 
tively [T7I. Our lattice QCD calculation, however, is 
missing some ingredients from the real world which 
means that we must adjust the experimental values we 
use in our calibration. The key missing ingredients are 
electromagnetism, c quarks in the sea and the possibility 
for the rjb and rjc mesons to annihilate to gluons, which 
we do not allow for in determining our yy^ and rjb corre- 
lators. These effects all act in the same direction, that 
of lowering the meson mass in the real world compared 
to that in our lattice QCD world. We estimate the total 
shift from these effects for the rjc to be -5.4(2.7) MeV and 
for the r]b, as -9(6) MeV The appropriate 'experi- 

mental' masses for the r]c and rji, for our calculations are 
then 2.986(3) GeV and 9.400(7) GeV. 

Since we need to allow for the three ingredients missing 
from our lattice QCD calculation when we determine me- 



son masses in section III we give in Table |II] a summary 
of our estimates of these effects jlOl [1^ . These estimates 
will be used to shift the lattice QCD results for compar- 
ison to experiment. Effects on decay constants are much 
smaller and we do not apply shifts in that case but simply 
include an additional uncertainty in the error budget. 

For consistency we fit a similar functional form to 
all quantities. This form must take account of physical 
heavy quark mass-dependence; discretisation errors and, 
for heavy-charm and heavy-strange mesons, mistuning of 
c and s quark masses. We use the standard constrained 
fitting techniques that we earlier applied to the correla- 
tors ^15^ . For the dependence on heavy quark mass and 



fiM,a) = fo{M) X 

7 3 1 
i=0 j,k=Q 1=0 

+Sfs + 5fc 



Mo 
M 



,aA, 



TT TT TT 



21 



(5) 



The quantity that we use for the heavy quark mass, M, 
is given by M = M^^ . /o is a function giving the 'lead- 
ing power' behaviour expected for each quantity. This 
is either derived from HQET or potential model expec- 
tations and takes the general form A{M/Mq)^. For the 
decay constants fn, and fn, we multiply this by the ra- 
tio of as values at the b and the c raised to the power of 
— 2//3o = -2/9 for rt/ — 3. This is the expected prefactor 
from resumming leading logarithms in HQET [inj. For 
Us we take ay from lattice QCD [1 (g^. We take Mq 
to have the value 2 GeV so that the factor Mq/M is ap- 
proximately IGeV/mfc. We tabulate the different forms 
for /o in Table |III| along with the prior values taken for 
A and b. b is allowed to float for the fit to /,,^ . In other 
cases it is fixed to the expected value but we have checked 
that allowing it to float returns the expected value within 
errors. We take the same prior for A of ± 2 in all cases. 

The sum to the right of the leading term in- 
cludes higher order corrections to the physical mass- 
dependence. These take the form of powers of Mq/AI, 
again using Mq — 2 GeV. We allow for 8 terms in the 
sum so that there is enough leeway to describe (by Tay- 
lor's theorem) any physically reasonable functional form 
in the fixed mass range from c to b. For the heavy-charm 
case we in fact fit from M — 4 GeV upwards so that 
the functional form is that appropriate to the unequal 
valence mass case. 

The other terms in the sum of eq. ([5| allow for sys- 
tematic errors resulting from sensitivity to the lattice 
spacing. Such discretisation errors depend on the lattice 
momentum cut-off, 7r/a, but can have a scale set by the 
different masses involved in the quantity under study. We 
allow for discretisation errors appearing with a scale of 
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mi and m2, where mi and m2 are the two quark masses 
in the meson (they will be the same in heavyonium). To 
be conservative we allow in addition further discretisa- 
tion errors with a scale of Aqcd where we take Aqcd — 
0.5 GeV. The powers of lattice spacing that appear in 
the terms must be even since discretisation errors only 
appear as even powers for staggered quarks. For the de- 
cay constants the Cijki are normalised so that cqooo = 1- 
For the mass differences the fits are normalised so that 
^ is 1 and cqooo floats. This is simply so that the fit 
can allow for significant discretisation errors when the 
physical mass difference is very small (particularly for 
the case of the Be to be discussed in section III C I . The 
prior values for the other Cijki are taken to be the same 
for all i, j, k and I but vary depending on the size of 
discretisation errors for the quantity being fit. They are 
larger for heavyonium than for heavy-strange quantities, 
for example. The values used are tabulated in Table [Till 
The mistuning of the strange and charm quark masses, 
where relevant, can be handled very simply because our 
tuning of these masses is in fact very good. We simply 
include an additional additive factor in the fit of 

ds 



Sfs 



M 

(™773,latt ~ ™77a,contnm) 

for heavy-strange mesons and 
dr 



^ami^^ ^ ^am2y 



(6) 



Sfc = (Cc 



^ ami ^ -\- 0-1^2 ^ ' 



(m^^ [att m^y^ contnm) 



(7) 



for heavy-charm mesons. The forms above allow for lin- 
ear quark mass dependence away from the tuned point. 
We do not need to include higher order terms because 
we are so close to the tuned point but we do allow for an 
M-dependent slope with discretisation errors (although 
in most cases neither of these additions makes any differ- 
ence) . 

To tune the strange quark mass we use the r^s, an un- 
physical ss pseudoscalar meson whose valence quarks are 
not allowed to annihilate. Lattice QCD simulations show 
that its mass m^^^contnm = 0.6858(40)GeV [TJ when the 
strange quark mass is tuned (from the K meson). Be- 
ing a light pseudoscalar meson, the square of its mass 
is proportional to the quark mass. To tune the c quark 
mass we use the r]c meson, as discussed earlier. The rjc 
meson is far from the light quark limit and so the meson 
mass is simply proportional to the quark mass. Cg and Cc 
are dimensionful coefhcients that represent physical light 
quark mass dependence and can be compared between 
lattice QCD calculations and with models. 

We do not include correlations between the results for 
different M on a given ensemble. We have not measured 
these correlations and the empirical Bayes criterion sug- 
gests that they are small. If we include a correlation 
matrix by hand for the results it makes very little dif- 
ference, a fraction of a standard deviation, to the final 
results. 
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1 0.66 


-0.244 


1.92020(16) 0.3044(4) 


0.81 


-0.335 


2.19381(16) 0.3491(5) 


0.825 


-0.344 


2.22013(15) 0.3539(5) 


0.85 


-0.359 


2.26352(15) 0.3622(5) 



0.44 

0.63 

0.66 

0.72 

0.753 

0.85 



-0.12 

-0.226 

-0.244 

-0.28 

-0.3 

-0.36 



1.42402(13) 
1.80849(11) 
1.86666(10) 
1.98109(10) 
2.04293(10) 
2.21935(10) 



0.21786(21) 
0.25998(20) 
0.26721(20) 
0.28228(22) 
0.29114(24) 
0.31900(27) 



0.3 

0.413 

0.43 

0.44 

0.45 

0.7 

0.85 



-0.06 

-0.107 

-0.115 

-0.12 

-0.125 

-0.27 

-0.36 



1.03141(8) 
1.28057(7) 
1.31691(7) 
1.33816(7) 
1.35934(7) 
1.86536(5) 
2.14981(5) 



0.15205(11) 
0.17217(11) 
0.17508(11) 
0.17678(11) 
0.17850(11) 
0.22339(12) 
0.25658(12) 



0.273 

0.28 

0.564 

0.705 

0.76 

0.85 



-0.0487 

-0.051 

-0.187 

-0.271 

-0.305 

-0.359 



0.89935(10) 

0.91543(8) 

1.52542(6) 

1.80845(6) 

1.91567(6) 

2.08753(6) 



0.11864(24) 
0.11986(21) 
0.16004(16) 
0.18071(16) 
0.18962(17) 
0.20576(16) 



0.193 

0.195 

0.4 

0.5 

0.7 

0.85 



-0.0247 

-0.02525 

-0.101 

-0.151 

-0.268 

-0.359 



0.66628(13) 

0.67117(6) 

1.13276(7) 

1.34477(8) 

1.75189(7) 

2.04296(7) 



0.0882(3) 

0.08846(11) 

0.1149(4) 

0.1260(5) 

0.1498(5) 

0.1708(6) 



TABLE IV: Results for the masses and decay constants in 
lattice units of the goldstone pseudoscalars made from va- 
lence HISQ heavy quarks on the different MILC ensembles, 
enumerated in Table |l] Columns 2 and 3 give the correspond- 
ing bare heavy quark mass and the e parameter, calculated 
at tree-level in m^a |TD]. This corresponds to a coefficient 
for the Naik 3-link discretisation correction of 1 + e. Meson 
masses from fitting these correlators using a simpler fitting 
form are given in ^ . Results given here are in agreement but 
somewhat more accurate. The results for heavy quark masses 
close to charm are also given in [10] . 



We also do not include effects from sea quark mass- 
dependence but, based on earlier work [10], we are able 
to estimate an uncertainty for that in our final results. 



III. RESULTS 



The correlators for pseudoscalar heavyonium mesons 
have very little noise and we can readily obtain ground- 
state masses with statistical errors in the fourth or fifth 
decimal place and ground-state decay constants with er- 
rors of 0.1%. Our results on each ensemble are given in 
Table HVl 



Results for /^^ are plotted against M^^ in Figure 
Discretisation errors are apparent in this plot and leai 
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FIG. 1: Results for the pseudoscalar heavyonium decay con- 
stant plotted as a function of the pseudoscalar heavyonium 
mass. Results for very coarse, coarse, fine, superfine and ul- 
trafine lattices appear from left to right. The colored dashed 
lines give the fitted function for that lattice spacing. The 
black line with grey error band gives the physical curve de- 
rived from our fit. The black circles with error bars at Af,,^ 
and Mr,^ are the values for the heavyonium vector decay con- 
stant at these physical points derived from the experimental 
leptonic widths for the J/i/) and T. The left-most black circle 
corresponds to the fictitious pseudoscalar 77^ particle whose 
decay constant was determined in [14) . 



to results at each value of the lattice spacing deviating 
substantially from the physical curve as the quark mass 
is increased. We fit the results to a physical curve al- 
lowing for discretisation errors as a function of the mass, 
as described in section |ll] and using the priors from Ta- 
ble III The power, 6, in eq. ([s]) is allowed to float in the 
fit. 

The obvious approach from which to gain some phys- 
ical insight in this case is that of the nonrelativistic po- 
tential model. In its simplest form this involves solv- 
ing Schrodinger's equation for the wavefunction of a two- 
particle system with reduced mass /it (= mb/2 for two h 
quarks) in a potential V(r) which is a function of the 
radial separation, r. At short distances we expect a 
Coulomb- like potential from QCD, and at large distances 
a string-like linear potential. However, other phenomeno- 
logical forms that interpolate between these two at in- 
termediate distances also work well at reproducing the 
bound-state spectrum, see for example |21j . The wave- 
function is useful for a first approximation in calculations 
of transition rates. In this sense, the wavefunction at the 
origin, "0(0)1 can be related to the decay constant by 
-0(0) = /^^^M,„712. However, 0(0) must be renor- 
malised before it can be related to a physical matrix 
element and some of the radiative corrections are very 
substantial [H]. In addition values of 0(0) vary widely 
with different forms for the potential that reproduce the 
same bound state spectrum because the spectrum itself 
provides little constraint on the potential at short dis- 
tances |22j . Here we will make comparisons of our lattice 



QCD results to those from potential models but it is im- 
portant to realise that the lattice QCD results for decay 
constants represent well-defined matrix elements in QCD 
and not model calculations. 

For a potential model with potential power count- 

oc ^3/(4+2Af) ^ggg^ gxam- 



ing arguments yield 0(0) 



pie, [13]). Then we would expect our fit for /^^ to need 
& = 1 for = — 1 but 6 = for iV = 1, the two extremes 
of the QCD heavy quark potential. Simply from com- 
paring values at c and h we might infer 6 « 0.5. In fact 
our fit gives the result h — —0.08(10) but with significant 
power corrections in 1/M, so that a simple power in M 
does not describe the results using our parameterisation. 
The physical curve that we extract of dependence on the 
heavyonium meson mass is shown as the grey band in 
Figure [T] 

The fit has of 1.2 for 29 degrees of freedom and 
allows us to extract results for c and h quarks. The re- 
sult for /^^ agrees within Icr with our earlier result of 
0.3947(22) "^GeV [lU] where we fit results at c only but 
included additional ensembles at different values of the 
sea ujd quark masses. Results for h quarks give: 



frib I f vc 



0.667(6) (2)GeV 
1.698(13)(5). 



(8) 



The first error comes from the fit and the second from ad- 
ditional systematic errors from effects not included in our 
lattice QCD calculation, i.e. electromagnetism, c quarks 
in the sea and (since we have not extrapolated to phys- 
ical u/d sea quark masses here) sea quark mass effects. 
Both errors are split into their component parts in the 
error budget of Table [V] We estimated the effects of elec- 
tromagnetism on fri^ from a potential model in TIT . We 
take the same 0.4% error for /^^^ since it is a more tightly 
bound particle but with smaller electromagnetic charges. 
There is then some cancellation of the effect in the ratio 
fvb/fvc- The effects of c quarks in the sea were shown 
to be similar to that of the hyperfine potential in [TU] 
and the effect on /^^ can then be estimated from the 
difference between /^^ and its associated vector particle. 
This is very small as we show below. We therefore ex- 
pect that missing c in the sea has a negligible effect on 
f,f^ and we estimate 0.2% on /^^ where it is magnified by 
(mil / rric)'^ . Sea quark mass effects on /^^ were shown to 
be very small in [10 , at the same level as the statistical 
errors of 0.1%. For /^^ we expect even smaller effects 
because it is a smaller particle. We take a 0.1% error 
nevertheless, but allow for some cancellation in the ratio 

of fvblfri.- 

The two rightmost black points (at and M^^) in 
Fig. [1] give the experimental values for the decay con- 
stants of the corresponding vector heavyonium mesons, 
J/0 and T, for comparison to the results calculated here 
in lattice QCD for the rjc and rji,. The decay constant for 
a vector meson can be defined by: 



< > /3 = fvmv 



(9) 
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Error 


fvb 




statistics 


0.6 


0.6 


M extrapoln 


0.2 


0.1 


a? extrapoln 


0.5 


0.4 


r\ 


0.4 


0.1 


r\/a 


0.5 


0.3 


Mr,, 


0.00 


0.05 


sea quark mass effects 


0.1 


0.05 


electromagnetism 


0.4 


0.2 


c in the sea 


0.2 


0.2 


Total (%) 


1.0 


0.9 



TABLE V: Full error budget for fr,^ and the ratio ft^^/frj, in 
%. See text for a fuller description of each error. The total 
error is obtained by adding the individual errors in quadra- 
ture. 



It has the advantage here that it can be extracted very 
accurately from experiment because vector heavyonium 
mesons can annihilate, through the vector currrent, to 
a photon, seen as two leptons in the final state. The 
relationship between the leptonic decay width and the 
decay constant is: 



T{Vh ^ e+e-) 



47r 



^QED^h 



my 



(10) 



where et is the electric charge of the heavy quark in units 
of e. The experimental results |17j give fj/^ = 407(5) 
MeV and /x — 689(5) GeV, remembering that the elec- 
tromagnetic coupling constant runs with scale and using 
l/aQED{m,) = 134 and l/aQED{mb) = 132[24. Thus 
1% accurate results for this decay constant are available 
from experiment, and can be used to test lattice QCD. 
Lattice QCD calculations of the T decay constant can be 
done |25j but they are not yet as accurate as the results 
we give here for the 7]^. 

The surprising result that we find on comparing the 
vector decay constant from experiment to the pseu- 
doscalar decay constant from lattice QCD is how close 
they are. In the nonrelativistic limit, where spin effects 
disappear, the vector and pscudoscalar become the same 
particle. Away from this point, however, there can be 
substantial relativistic corrections, particularly for char- 
monium. Instead we find that the pseudoscalar decay 
constant is 3% lower than the vector in both cases with 
an error of 1-2%. 

Unfortunately this cannot be directly tested through 
decay modes of the r]c or ry^. The decay rate to two 
photons is indirectly related to the decay constant as the 
leading term in a nonrelativistic approximation: 



(11) 



This formula has radiative and relativistic corrections at 
the next order. The decay width is not known for the rji, 
and only very poorly known for the rj^, with the Particle 
Data Group estimate given as 7.2(2.1) keV [T^. Substi- 



tuting this into eq. (11) and taking rric = justi- 
fiable at this order, gives /r,^ = 0.4(1) GeV, where only 
the large error from experiment is shown. This is consis- 
tent with our value but much less accurate so does not 
provide a useful test. 

As discussed earlier, a direct comparison of lattice 
QCD results for /^^ and potential model values for ip{0) 
is not particularly useful. Values for V'(O) for the ground 
state in bottomonium vary by a factor of 1.5 for different 
forms for the potential in (55]. This variation is reduced 
somewhat, and radiative corrections cancel, if we com- 
pare the ratio of values at b and c. Here the lattice QCD 
result above of 1.698(14) favours the strong variation of 
-0(0) with quark mass seen in the Cornell potential. For 
this potential [22* gives a ratio V'()(0)/V'c(0) of 3.1, yield- 
ing a decay constant ratio of 1.8. 

Figure [l] also includes as the leftmost black point a 
value for the decay constant of the rjs as determined from 
lattice QCD [TJ. Although our fit becomes unstable be- 
low M of 2 GeV, it is interesting to see that /^^ does 
not look out of place on this plot as the light and heavy 
sectors are smoothly connected together. 



B. rnB^ and /s^ 

Our calculations for heavy-strange mesons were de- 
scribed in [3| and so we only add briefly to that discussion 
here. In Table [Vl] we give our full set of results, including 
values at a variety of strange quark masses for complete- 
ness. In [3] we used the heavy-strange mass itself as a 
proxy for the heavy quark mass and obtained good agree- 
ment for the mass of the Bg with experiment and a value 
for fs, of 225(4) MeV. 

Here, for consistency with the other calculations, we 
use instead M^^ for the heavy quark mass and the fit 
form given in eq. ([5|. For the heavy-strange meson mass, 
as in [3] , we fit to the mass difference: 



(12) 



We take account of mistuning of the strange quark mass 
using the factor given in eq. For the decay constant 
fit we fix the power of the leading M-dependence, b = 
—0.5. Allowing b to float gives results for b in agreement 
with this value to within 20%. 

Our fit to AH^,hh& is shown in Figure [2] and gives 
= 0.2 for 17 degrees of freedom. The values ex- 
tracted at the c and b masses agree well, within Icr, 
with our earlier results [ll[TUj. When account is taken of 
electromagnetic and other effects missing in the lattice 
calculation these earlier results translate into values for 
mo, = 1.969(3) GeV [ID] and ms, = 5.358(12) GeV g]. 
The increased error at the b results from increased sta- 
tistical and discretisation errors for heavier quark masses 
as well as the extrapolation in M. Our result for 
agrees within the 12 MeV error with that deter- 
mined from full lattice QCD using a completely different 
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Set 


rusa Mrj^a mha 


MhsO. 


fa, a 




1 


0.061 0.50490(36) 0.66 


1.3108(6) 


0.1913 


;7) 




0.81 


1.4665(8) 


0.1970 


;io) 




0.066 0.52524(36) 0.66 


1.3164(5) 


0.1929 


'J) 




0.825 1.4869(7) 


0.1994 


;io) 


2 


0.0492 0.41436(23) 0.44 


0.9850(4) 


0.1500 


;5) 




0.63 


1.2007(5) 


0.1559 


'J) 




0.85 


1.4289(8) 


0.1613 


;io) 




0.0546 0.43654(24) 0.44 


0.9915(4) 


0.1516 


;5) 




0.66 


1.2391(5) 


0.1586 


;6) 




0.85 


1.4348(7) 


0.1634 


;9) 


3 


0.0337 0.29413(12) 0.3 


0.70845(17) 


0.1054 


;2) 




0.413 


0.84721(23) 


0.1084 


'2) 




0.7 


1.1660(4) 


0.1112 


;5) 




0.85 


1.3190(5) 


0.1123 


;6) 




0.0358 0.30332(12) 0.3 


0.71119(16) 


0.1061 


'2) 




0.43 


0.86982(23) 


0.1094 


'2) 




0.44 


0.88152(23) 


0.1096 


'3) 




0.7 


1.1684(4) 


0.1121 


;4) 




0.85 


1.3214(5) 


0.1131 


;6) 




0.0366 0.30675(12) 0.3 


0.71223(16) 


0.1063 


;2) 




0.43 


0.87079(22) 


0.1097 


'2) 




0.44 


0.88249(23) 


0.1099 


;3) 




0.7 


1.1694(4) 


0.1124 


;4) 




0.85 


1.3223(5) 


0.1135 


;6) 


4 


0.0228 0.20621(19) 0.273 0.59350(24) 0.0750 


'3) 




0.564 0.9313(5) 


0.0754 


'6) 




0.705 1.0811(8) 


0.0747 


;8) 




0.85 


1.2279(10) 


0.0742 


;io) 


5 


0.0161 0.15278(28) 0.193 


0.43942(33) 


0.0553 


'4) 




0.5 


0.8027(10) 


0.0541 


;i2) 




0.7 


1.0152(18) 


0.0513 


;22) 




0.85 


1.1657(24) 


0.0495 


'30) 




0.0165 0.15484(14) 0.195 


0.44270(28) 


0.0555 


;3) 




0.5 


0.8038(8) 


0.0546 






0.7 


1.0169(12) 


0.0526 


;i6) 




0.85 


1.1684(16) 


0.0517 


;2i) 



TABLE VI: Results for the masses and decay constants in 
lattice units of the goldstone pseudoscalars made from valence 
HISQ heavy quarks with valence HISQ strange quarks on the 
different MILC ensembles, enumerated in Table [l] Column 2 
gives the s mass in lattice units, with several values on some 
ensembles around the correctly tuned value. Column 3 gives 
the corresponding mass for the goldstone pseudoscalar made 
from the s quarks, which is used for tuning. Column 4 gives 
the heavy quark mass. The corresponding values of the Naik 
coefficient are given in Table |IV| Many of these results were 
given earlier in j?) [TU] . 



method (NRQCD) for the b quark [18] with very differ- 
ent systematic errors, providing a stringent test of lattice 
QCD. Our results also agree well with experiment [T7] 
{iTiDs = 1.968 GeV and niB^ = 5.367 GeV) and this pro- 
vides a very strong test of QCD. 

The fit to the decay constant, fn^ , is shown in Figure[3] 
and gives — 0.3 for 17 degrees of freedom. Again 
results at the b and c agree within la with our earlier 
results which are: /z,, = 0.2480(25) GeV dU] and /s, = 
0.225(4) GeV [4]. 

Figures [2] and |3] give the physical fit curves as a func- 



0.70 




2 4 6 8 10 

(GeV) 



FIG. 2: Results for the difference, An^Mh between the heavy- 
strange pseudoscalar meson mass and one half of the pseu- 
doscalar heavyonium mass. Results for very coarse, coarse, 
fine, superfine and ultrafine lattices appear from left to right. 
The lattice QCD results have been adjusted for slight mistun- 
ing of the s quark mass. The colored dashed lines give the 
fitted function for that lattice spacing. The black dashed line 
with grey error band gives the physical curve derived from 
our fit. The black circles with error bars at Af^^ and Af^^ are 
the experimental values adjusted for the effects from electro- 
magnetism, rib/rjc annihilation and c quarks in the sea, none 
of which is included in the lattice QCD calculation. 
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FIG. 3: Results for the pseudoscalar heavy-strange decay con- 
stant plotted as a function of the pseudoscalar heavyonium 
mass. Results for very coarse, coarse, fine, superfine and ul- 
trafine lattices appear from left to right. The lattice QCD 
results have been adjusted for slight mistuning of the s quark 
mass. The colored dashed lines give the fitted function for 
that lattice spacing. The black dashed line with grey error 
band gives the physical curve derived from our fit. 



tion of Mjj^ . As expected, the curves are very similar 
to those in [1] since to a large extent the change is sim- 
ply a rescaling of the x-axis. However they provide a 
consistency check that the parameterisation we use here, 
taking a different quantity to represent the heavy quark 
mass, works just as well. 
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Set rrica 


ruha 


Mn^a 


fa, a 


2 0.63 


0.85 


2.01651(10) 0.2854(2) 


3 0.413 


0.7 


1.57733(7) 


0.1916(2) 




0.85 


1.72373(6) 


0.2004(1) 


0.43 


0.7 


1.59489(7) 


0.1938(2) 




0.85 


1.74105(6) 


0.2030(1) 


0.44 


0.7 


1.60522(6) 


0.1952(1) 




0.85 


1.75122(6) 


0.2044(1) 


4 0.273 


0.564 


1.21799(8) 


0.1329(2) 




0.705 


1.36350(8) 


0.1367(2) 




0.76 


1.41872(8) 


0.1380(2) 




0.85 


1.50727(8) 


0.1402(2) 


0.28 


0.564 


1.22562(8) 


0.1338(2) 




0.705 


1.37103(8) 


0.1376(2) 




0.76 


1.42621(9) 


0.1390(2) 




0.85 


1.51471(9) 


0.1413(2) 


5 0.195 


0.4 


0.90566(8) 


0.0967(3) 




0.5 


1.01457(9) 


0.0985(4) 




0.7 


1.22392(10) 0.1005(4) 




0.85 


1.37366(10) 0.1018(5) 



TABLE VII: Results for the masses and decay constants in 
lattice units of the goldstone pseudoscalars made from va- 
lence HISQ heavy quarks with valence HISQ charm quarks 
on the different MILC ensembles, enumerated in Table |T] Set 
1 is missing because rrica is already close to the highest heavy 
quark mass that we use. Column 2 gives the c mass in lat- 
tice units, with several values on some ensembles around the 
tuned c mass, and column 3 the heavy quark mass. The corre- 
sponding values of the Naik coefficient are given in Table |IV[ 

C. rriB^ and fs^ 

Heavy-charm mesons are of interest because a family of 
gold-plated be mesons exists of which only one, the pseu- 
doscalar Be [26t [27] , has been seen. Traditionally these 
particles have been viewed as further examples, beyond 
bb and cc, of a heavy-heavy system and therefore a test 
of our understanding of this area, be mesons, however, 
have a lot in common with heavy-light systems. In fact 
they provide a bridge between heavy-heavy and heavy- 
light and so test our contol of QCD much more widely. 
The more accurately we can do these tests, the better 
they are. 

Lattice QCD calculations of the Be mass can be done 
very accurately. Indeed the mass of the B,. was pre- 
dicted ahead of experiment with a 22 MeV error [2H] 
using NRQCD for the b quark and the 'Fermilab' clover 
action for the c quark. The error was later reduced to 10 
MeV by using a more highly improved action, HISQ, for 
the charm quark [TS]. Here we use the HISQ action for 
both the c quark and the heavier quark up to the b mass 
to obtain results in a completely different heavy quark 
formalism. In addition we calculate the decay constant 
of the Be for the first time in full lattice QCD. 

To determine the Be mass we use the mass difference 
to the average of the associated heavyonium states: 

AH^,hh = Mh^ - ^(M,^ + Af,J. (13) 




Mr^, (GeV) 



FIG. 4: Results for the mass difference between the He meson 
and the average of the associated heavyonium pseudoscalar 
meson masses plotted as a function of the pseudoscalar heavy- 
onium mass. Results for coarse, fine, superfine and ultrafine 
lattices appear from left to right. The lattice QCD results 
have been adjusted for slight mistuning of the c quark mass. 
The colored dashed lines give the fitted function for that lat- 
tice spacing. The black line with grey error band gives the 
physical curve derived from our fit. The black circle with 
error bar at M^^^ gives the experimental value adjusted for 
the effects from electromagnetism, r]b/r]c annihilation and c 
quarks in the sea, none of which is included in the lattice 
QCD calculation. 
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FIG. 5: Results for the mass difference between the He meson 
and the average of the associated heavyonium pseudoscalar 
meson masses plotted as a function of the pseudoscalar char- 
monium mass. Results are given for two heavy quark masses 
on fine lattice set 3 (pink bursts) and four heavy quark masses 
on superfine lattices set 4 (green crosses) . Lines are drawn to 
guide the eye. 



^Ha,hh is a measure of the difference in binding energy 
between the symmetric heavyonium states made of c and 
h quarks and the heavyonium state made of two different 
mass quarks, c and h. Here we map out AH^,hh as a func- 
tion of the heavy quark mass, and reconstruct Mb^ from 
'-^H^,hh determined at h = b. An^^hh can be determined 
with high statistical accuracy because all of the states 
involved have very little noise. The fact that Ah^^hh is 



10 



very small ( for = rric by definition and less than 
100 MeV when ruh = mt,) also means that lattice errors 
from, for example, the uncertainty in the lattice spacing 
are very small. In fact, for this calculation, as discussed 
below, key sources of error are the uncertainties from 
electromagnetic, annihilation and c-in-the-sea shifts to 
the masses. 

Table |VII| gives our results for the masses and decay 
constants of the He mesons calculated using quark masses 
that are close to that of the c quark mass on each ensem- 
ble and then all the heavier masses for h. We give re- 
sults for more than one value of the c quark mass on the 
fine and superfine ensembles (sets 3 and 4) so that slight 
mis-tuning in the c quark mass can be corrected for. It 
is clear from the results that l^^u^.hh can be calculated 
with a statistical accuracy of better than 1 MeV. Errors 
from uncertainties in the lattice spacing are also at this 
level. 

Figure H] shows An^^hh plotted against A/^^ for the re- 
sults at different values of the lattice spacing. A fairly 
clear linear dependence is evident. An^^hh would be ex- 
pected to increase linearly with M^^ at large Af^^ , in the 
same way as An^.hh, from a simple potential model argu- 
ment. The binding energy of the rjh becomes increasingly 
negative, roughly in proportion to Af^,^ as it increases (at 
least for a potential with N = —1), whilst the binding 
energy of the He meson does not change. A corollary of 
this is that the dependence of An^^hh on (as proxy 
for TTic) would also then be expected to be linear with a 
slope of opposite sign and roughly three times the magni- 
tude. The factor of three is because the binding energy of 
the ric becomes more negative as Af^^ increases, with the 
same dependence as the rih binding energy has on Af^^. 
The He binding energy will also become more negative 
but have double the slope because the reduced mass of 
the He system is roughly rric rather than nie/1 for the 
rjc- On top of this Af,,^ appears halved in An^^hh- 

Interestingly this factor of -3 does seem to be approxi- 
mately true in comparing Figure [5j which shows the de- 
pendence of An^.hh on M^^, with Figurelil FigureHlgives 
a slope of « 0.012 (over the full range) and Figure fsl gives 
slopes varying from -0.03 to -0.04 over a small range in 
Af^^, as Mrj^ increases. In our fit to An^.hh we include 
the effect of mistuning rric (from eq. Q) and obtain con- 
sistent values from that. 

We fit An^^hh as a function of A/^,^ (above 4 GeV) us- 
ing the fit form described in section [ll] The leading mass 
dependence is taken to be Af^^ — Af^^, so that An^^hh 
vanishes when Af^^ = M^^ as it must by definition. As 
described in section [IT] we include a sum of power cor- 
rection terms and lat tice spacing dependent terms with 
priors given in Table III The fit gives of 0.3 for 11 



degrees of freedom and result: 

Aij,,b6 = 0.065(9)GeV. 



(14) 



The resulting physical curve of heavy quark mass de- 
pendence is shown in grey on Figure |4] The comparison 
to experiment is given by the black dot with error bar at 



Error 






ABc,bs 


statistics 


8.4 


0.7 


0.5 


M extrapoln 


3.1 


0.2 


0.2 


a? extrapoln 


10.9 


0.7 


0.4 


ri 


0.7 


0.6 


0.3 


ri/a 


1.4 


0.8 


0.3 




0.9 


0.5 


0.3 


sea quark mass effects 


1.5 


0.1 


0.1 


electromagnetism 


3.1* 


0.4 


0.2* 


c in the sea 


5.3* 


0.04 


0.1* 


rib,c annihiln 


2.7* 






Total (%) 


18 


1.6 


0.9 



TABLE VIIL Full error budget for AB,,bb, /s, and Ag^bs 
given as a percentage of the value. See the text for a fuller de- 
scription of each error. The total error is obtained by adding 
the individual errors in quadrature, except for the final three 
systematic errors (starred) for Ag^^bb and AB^,bs which are 
correlated and so simply added together before being com- 
bined in quadrature with the others. 



h = b. This experimental result has been shifted to be 
the appropriate value to compare to our lattice QCD cal- 
culation as we now describe. The current world-average 
experimental result for Mb, - 0.5(A/^^ -I- Af^J is 92(6) 
MeV [17] . There is a sizeable experimental error coming 
mainly from the Be but also from the ryt,. Our lattice 
QCD calculation is done in a world without electromag- 
netism or c quarks in the sea and in which the rji, and 
rjc do not annihilate. The absence of these effects (i.e. 
to compare to our lattice result) produces shifts to the 
masses as discussed in section |lll Estimated values for 
the shifts are given in Table [TTj The net effect is to shift 
the experimental value of As^^hh down by -8(7) MeV, 
where the error takes the shifts to be correlated. The 
'experimental' value of AB^.hh to compare to our lattice 
result is then 84(9) MeV, marked on Figure |4j Our lat- 
tice result agrees with experiment, once these shifts are 
made, within 2a. 

From As^fib we can reconstruct the Be meson mass, 
now applying the shifts above to the lattice QCD calcu- 
lation to obtain a result that can be compared to exper- 
iment. This gives the result 



Mb, = 6.259(9) (7)GeV. 



(15) 



Here the first error comes from the fit and the second er- 
ror from the shifts applied to include missing real world 
effects as well as experimental uncertainties in the rj^ and 
r}c masses. As can be seen, this is a sizeable part of the 
total error in this case. We also include in this second 
error an estimate of sea quark mass effects using results 
from |10] . There we saw no such for an equivalent quan- 
tity for rriB)^ within 1 MeV statistical errors and so take 
that as the error here. Table |VIII| gives the complete 
error budget for As^fib breaking down both errors into 
their components. 

Our result for Mb, can be compared to experiment 
(6.277(6) GeV) and to our result from lattice QCD us- 
ing a completely different formalism, NRQCD, for the b 
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FIG. 6: Results for the heavy-charm decay constant plot- 
ted as a function of the c quark mass, given by the mass of 
the r]c meson. Results are given for multiple heavy quark 
masses on fine lattices (pink bursts) and superfine lattices 
(green crosses). Lines are drawn to guide the eye. 



quark (6.280(10) GeV [H]). We agree, within 2a with 
both results even allowing for the fact that the compar- 
ison within lattice QGD can be done before any shifts 
are made or errors allowed for them. This is a strong 
confirmation of the control over errors that we now have 
in lattice QGD. 

The method given here for determining tob^ (as for 
the method for tub^ 



in section III B ) does depend on the 
experimental rjb mass; the mass difference determined in 
lattice QGD is not particularly sensitive to it but when 
the mass is reconstructed from the difference, is 
added in. Recent results from the Belle collaboration [25] 
have = 9.402(2) GeV, significantly higher than the 
previous world-average [T7]. Using the Belle result for 
Mjj^ pushes our values for niB^ and m^^ 6 MeV higher. In 
both cases this improves the agreement with experiment 
but is not significant given the 11 MeV error. Note that 
our earlier NRQGD results are hardly affected at all by 
a change in the rji, mass because they determined a mass 
difference to the spin-average of the T and rji, masses, 
which is dominated by the T mass. 

Results for the He decay constant, fn^, are also given 
in Table [VlTl The rate for Be leptonic decay to liy via a W 
boson is proportional to the square of the decay constant 
multiplied by GKM element Vcb as in eq. Q. In prac- 
tice this decay will be very hard to see experimentally, 
but a lattice QGD calculation of the decay constant also 
provides a useful test for phenomenological model calcu- 
lations. 

The results at different values of ttLc can again be used 
to tune the decay constant accurately to the result at the 
physical c quark mass. Figure [6] shows the dependence 
of /ifc on Mjj^ acting as a proxy for the c quark mass. 
Results on fine and superfine lattices are shown - there 
is clear agreement on the physical slope of fn^ with M^^ 
between superfine and fine and it does not vary with the 
heavy quark mass. The slope is small, approximately 
0.06, but clearly visible. We will compare this to the 



FIG. 7: Results for the heavy-charm decay constant plotted 
as a function of the pseudoscalar heavyonium mass. Results 
for coarse, fine, superfine and ultrafine lattices appear from 
left to right. The lattice QGD results have been adjusted for 
slight mistuning of the c quark mass. The colored dashed lines 
give the fitted function for that lattice spacing. The black line 
with grey error band gives the physical curve derived from our 
fit. 



slope for fn^ with rris in section IV 

The H(. decay constant is plotted as a function of M^^ 
in Figure [7] Notice that it is much flatter than the cor- 
responding plot for /^^ (Figure [T]). We expect behaviour 
as 1/ y/Mrjf^ whether we view heavy-charm as a heavy- 
light system (in which case the behaviour will be similar 
to heavy-strange) or as a heavy-heavy system (in which 
case the argument becomes that tp{0) depends on the 
reduced mass /i, tending to nic for large nih, and then 
the decay constant falls as the square root of the heavy 
mass) . 

As before, we fit fn^ to the function of Af^^ (above 4 
GeV) described in section [llj We take the leading term 
given in Table |III| to be that expected from HQET ar- 
guments appropriate to heavy-light physics. Our fit has 
— 0.7 for 11 degrees of freedom and gives result: 



/s^ - 0.427(6) (2)GeV. 



(16) 



Here the first error is from the fit and the second from ad- 
ditional systematic effects missing from our lattice QGD 
calculation. These we estimate based on the arguments 
given for the rih in section |III A| The error from miss- 
ing electromagnetism and from sea quark mass effects 
we take to be the same as for the rjb at 0.4% and 0.1% 
respectively; missing c in the sea should be a factor of 
mc/rrib smaller at 0.04%. Table VIII gives the complete 
error budget. 

can be converted into a branching fraction for lep- 
tonic decay using the formula of eq. Q and the unitarity 
value of Vcb- We predict a branching fraction to tv of 
0.0194(18). The error here comes mainly from the exper- 
imental determination of the Be lifetime with a smaller 
effect from the uncertainty in Vcb- Our value for Jb^ con- 
tributes a 3% error. Because of helicity suppression the 
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branching fraction smaller for other lepton final states 
(8 X 10~^ to ^ly, for example). 

The value we obtain for fg^ can be compared to results 
from potential models. As discussed earlier in the context 
of , potential model results have a lot of variability and 
raw values for ip{0) need renormalisation. A more useful 
comparison is to compare ratios. Our lattice QCD results 
giye JbJU = 1-08(1) ^^df^JfB^ = 1.57(2). The range 
of potentials considered in [32] give values from 0.90 to 
1.02 for JbJU and 1.34 to 1.72 for f^Jfe^. Again 
the largest number is always from the Cornell potential. 
Potential model values for "0(0) converted to Jb^ simply 
using / = V(0)\/12/Ms^ yield results varying from 0.5 
to 0.7 GeV i.e. significantly larger than the well-defined 
value for /^^ from lattice QCD. 

The values for /^^ from potential models provide input 
to estimates of the production cross-section of the at 
the LHC. In the factorisation approach the cross-section 
is proportional to the square of /s^, with typical values 
for Jb^ being taken as 0.48 GeV [30j. Our results indicate 
that this could be leading to a 25% overestimate of the 
production rate. 



IV. DISCUSSION 

An interesting issue is to what extent the Be meson 
is a heavy-heavy particle and to what extent, a heavy- 
light one at the physical values we have for b and c quark 
masses. Here we address this by comparing the behaviour 
of Be properties to those of rjf, and Bg using the results 
from section ITlll 

An alternative to calculating An^^hh to study the 
heavy-charm meson mass is to take differences be- 
tween heavy-charm and heavy-strange and charm- 
strange mesons. We define 



A 



H^,hs 



Mr, -M 



(17) 



SO that ABc,hs is a positive quantity. Once again it 
amounts to a difference in binding energies but now be- 
tween a set of mesons that are all effectively 'heavy-light' 
states. Indeed a study of An^^hs shows us to what extent 
the Be can be considered a heavy-light particle rather 
than, or as well as, a heavy-heavy one. 

Figure [8] shows Ah^^hs, with all results tuned accu- 
rately to the correct c and s masses, as a function of the 
heavy quark mass, again given by the rjh mass. In fact 
AH^,hs shows very little dependence on the heavy quark 
mass above a value of M^^ of about 6 GeV. HQET would 
expect the leading m/i-dependent piece of AB^,hs to be 
given by the difference of the expectation values of the 
kinetic energy operator, p^/2rn/i, for the heavy quark in 
a heavy-charm meson and a heavy-strange meson, ignor- 
ing the effect of spin-dependent terms which are expected 
to be smaller. Figure [8] shows that this difference is not 
large i.e. the charm quark is behaving in a similar way to 
a light quark (but does have a larger expectation value for 
its kinetic energy operator as might be expected) when 



combined with a heavy quark of order twice its mass or 
heavier. 

We fit TTiHchs as described in section [ll| and using the 
fit form and priors tabulated in Table Illip Our fit has 
of 0.3 for 14 degrees of freedom. It returns the coefficient 
of the first term in M,J^^^ as -0.4(8)GeV/Af,,^ . This quan- 
tifies the statement made above about the slope of 1/ rrih 
corrections. The coefficient is not very accurately de- 
termined because we allow for many higher order terms. 
In fact the sign of the slope is clear from Figure [8] with 
a positive slope with M^^ corresponding to a negative 
value for the coefficient of the 1/M term, as expected. 

The variation of A^^^hs with M^^ agrees well with that 
found in our calculation using NRQCD b quarks [TS] giv- 
ing a slope of 0.07 at the b. Likewise the variation with 
Af,^^ also agrees well with the slope of 0.4 found in [18,. 

Our fit to AHc, hs is independent of our earlier fit to 
^H^.hh (although it uses some of the same numbers) and 
so the results provide a consistency check. We find at 
h = b that: 



A 



1.052(9)(3)GeV 



(18) 



which agrees well within la with the same quantity cal- 
culated using NRQCD b quarks [TH]. The result when 
h = cis consistent within la with double the result from 
rriB)^ — m,,^/2 given in [10,. The first error above is from 
the fit and the second from the systematic error for sea 
quark mass effects, taking the same 1 MeV as for As^,6f,, 
and the effects of missing electromagnet ism and c in the 
sea. The shifts and errors for these latter effects are given 
m Table HD and we take those errors to be correlated. 
The value above for AB^^bs combined with experimental 
results for Mb, and Md/FTT! (the net shift from Table |ll] 
amounts to a negligible 0.2 MeV) gives: 



Mb, = 6.285(9)(3)GeV, 



(19) 



consistent within 2a with our result from As^^bb given 
in section [nij and slightly more accurate. We therefore 
adopt it as our final result here. The complete error 
budget for AB,,bs is given in Table VIII 

In figure [9] we show the ratio of rjh and He decay con- 
stants to that of the , plotted from our physical curves 
as a function of M^^ . The ratio /^^ / fn, rises strongly 
with Mri^ , because of the big difference in the dynamics 
of heavy-heavy and heavy-strange mesons, whereas the 
ratio fHa/fHs tends to a constant at large M^^. As ex- 
plained in section |III C| this latter behaviour would be 
expected whether the heavy-charm is viewed as a heavy- 
heavy or heavy-light state, because the reduced mass of 
the heavy-charm system is controlled by the charm mass 
in the large heavy mass limit. 

Further insight comes from comparing the dependence 
of the heavy-charm and heavy-strange decay constants 
on rric and respectively. Figure [lO| plots the relative 
change of Jh, or Jh, to its value at the tuned mass point 
for a given relative change in the light quark mass. The 
strange quark mass is monitored by the value of M^^ , 
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FIG. 8: Results for the mass difTerence between the heavy- 
charm meson and the corresponding heavy-strange and 
charm-strange mesons plotted as a function of the pseu- 
doscalar heavyonium mass. Results for coarse, fine, superfine 
and ultrafine lattices appear from left to right. The lattice 
QCD results have been adjusted for slight mistuning of the c 
and s quark masses. The colored dashed lines give the fitted 
function for that lattice spacing. The black line with grey 
error band gives the physical curve derived from our fit. The 
black circles with error bars at M^^ and M,,^ give experi- 
mental values adjusted for the effects from electromagnetism, 
rjb/rjc annihilation and c quarks in the sea, none of which is 
included in the lattice QCD calculation. 
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FIG. 9: Results for the ratio of pseudoscalar decay constants, 
heavy-charm and heavy-heavy to heavy-strange plotted as a 
function of the pseudoscalar heavyonium mass. The results 
are obtained from the physical curves given in Figures [l] |3] 
and [3 



the charm mass by M^^. The results come from the fine 
lattices, set 3, where we have multiple rric and rus values 
close to the tuned point. Results are plotted for two 
values of the heavy quark mass, m/jO = 0.7 and ra^a — 
0.85 but little difference between them is seen. 

The dependence of Jh^ on is not very strong 
as expected since Jh, and fn differ only by around 20% 
for a change by a factor of 27 in light quark mass. The 
dependence of Jh^ on nic is larger by about a factor of 



FIG. 10: Comparison of the effect of 'detuning' the charm and 
strange quark masses on the heavy-charm and heavy-strange 
decay constants. Open squares show the fractional change in 
Jh^ for a given fractional change in M^^ (as proxy for rric) 
for two different heavy quark masses (in blue ruha — 0.7 and 
red mftfl = 0.85) on the fine lattices set 3. Burst show the 
fractional change in for a given fractional change in M^^ 
(as proxy for nis) for the same two heavy quark masses on set 
3. Lines are drawn to guide the eye. 



fvh 



fn 



A, 



A, 



0.394(2) 
0.452(2) 
0.501(3) 
0.546(4) 
0.586(4) 
0.623(5) 
0.655(6) 



0.249(2) 
0.251(2) 
0.249(2) 
0.244(3) 
0.237(3) 
0.231(4) 
0.224(4) 



0.417(6) 
0.427(3) 
0.434(4) 
0.435(4) 
0.433(5) 
0.429(6) 



0.477(2) 
0.520(3) 
0.554(4) 
0.581(6) 
0.605(7) 
0.626(9) 
0.645(11) 



0.000(0) 
0.004(1) 
0.015(1) 
0.027(2) 
0.039(3) 
0.050(5) 
0.061(8) 



0.956(6) 
0.994(7) 
1.014(6) 
1.028(6) 
1.038(7) 
1.045(8) 
1.050(8) 



TABLE IX: Values for the various quantities that we fit here 
evaluated at masses, M^^ , between that of c and b. These are 
obtained from our fit functions at a = and tuned s and c 
masses. All numbers are in GeV. There is no result for fn^ 
at 3 GeV because that point is not included in that fit. 



two. However the slope of the Figure 10 is 1/3 (see also 
Figure |6]), much less than the slope of 1 expected if /b^ oc 
nic- This latter behaviour would be approximately that 
expected in a heavy-heavy picture in which oc ^, 

with the reduced mass, fi, close to rric in the Be case. 
The linear behaviour of V'(O) would be consistent with the 



picture we have of the rjf, in Figure [T] where ^ ; 
using b « 0.5. 



M,,/4, 



V. CONCLUSIONS 

By using a relativistic approach to heavy quarks 
(HISQ) which has relatively small discretisation errors 
we have been able to map out the dependence on heavy 
quark mass of the pseudoscalar heavyonium, heavy- 
strange and heavy-charm decay constants and the heavy- 
strange and heavy-charm meson masses, complementing 
results in OH]. 

We find the heavyonium decay constant surprisingly 
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FIG. 11: Summary of heavy quark mass dependence of decay 
constants for the pseudoscalar He, Hs and rjh mesons. The 
grey bands show our physical curves from Figures [l] [S] and [7] 



close in value to the experimental results for the charnio- 
nium and bottomonium vector decay constants. Work 
is underway to confirm this result using NRQCD for the 
heavy quark and to establish accurate results for the cor- 
responding vector decay constants in lattice QCD. Al- 
though the r]ii decay constant has no simple connection 
to an observed experimental rate, it is useful for compar- 
ison and calibration of lattice QCD calculations in heavy 
quark physics since it can be determined to 1%, as we 
have done here. 

Our result for the Be meson mass agrees well using the 
two different mass splittings, AB^,hh and AB^,hs and also 
agrees with the experimental value. This is confirmation 
of our earlier result [TS] using NRQCD b quarks and HISQ 
light quarks. 

We determine the Be decay constant as 427(6) MeV, 
for the first time in full QCD, predicting a leptonic 
branching ratio for the Be to Tiy of 1.9(2)% (where the 
uncertainty comes from ts^, not fs^)- Our result for 
Jb^ is significantly smaller than that from some poten- 
tial model calculations, including those being used to es- 
timate LHC production cross-sections [3D] . The best way 
to determine the Be leptonic decay rate, and hence fs^, 
from experiment may be using a high luminosity e'^e" 
collider operating at the Z peak (SHI ■ 

By mapping out the dependence on the heavy quark 
mass of the He, Hg and ry/j decay constants we are able 
to see the differences between the three systems. This 
is summarised in Figure [TT] where we give the physical 
curves determined from our fits. In section [TV] we provide 
evidence that the Be behaves, at least in some ways, more 
like a heavy-light system than a heavy-heavy one. We 
previously noticed this effect in [38j when finding that 
the mass difference between B* and Be was very close to 
the difference between B* and Bg. 

Table |IX| gives results extracted from our fits at in- 



termediate values of from M^^ 



to M^i^ for com- 



parison to future lattice QCD calculations or to phe- 
nomenological models. The values are determined by 
evaluating our fit function in the continuum limit and 
at tuned s and c masses, corresponding to the black line 
in Figs.[lJ[2l|3l|4l[7]and[8l 

In Figure [T2| we summarise the current picture for the 
decay constants of gold-plated mesons, determined from 
lattice QCD and from experiment. For lattice QCD we 
use the best existing results which dominate the world 
averages [H H [lOl HH EB |32] . For the experimental val- 
ues for the unflavored vectors we use leptonic widths to 
e+e~ from the Particle Data Tables [T7j and eq. (10 1. 
For the flavored pseudoscalars the determination of the 
decay constant from experiment requires the input of a 
value for the associated CKM element, for example from 
the unitarity fit to the CKM matrix jTT] . We update the 
D and Dg experimental determinations to the averages 
including new results from BESIII [31] and Belle [3S] re- 
spectively. 

This plot goes beyond the traditional plot of the mass 
spectrum [T] to look at a number which is related to the 
internal structure of the meson. The energy scale for de- 
cay constants is controlled by internal momenta inside 
the meson and so is much compressed over the scale for 
masses (which covers a large range simply because quark 
masses have a large range). The pseudoscalar meson de- 
cay constants are well filled in but more work is needed 
to obtain the vector decay constants to the same level of 
accuracy. This is underway and once complete, this plot 
will provide a very stringent test of QCD that would be 
impossible with any method other than lattice QCD. 

From our results here and in [3j H] we see that the rela- 
tivistic heavy quark approach using the HISQ formalism 
can successfully give results for the h quark. Future work 
will use even finer lattices. For a = O.OSfm, for exam- 
ple, the b quark mass in lattice units is around 0.5 and 
so we can easily achieve this mass without the need for 
extrapolation. 
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